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We study the dynamics of entanglement in the XY spin chain subject to a staggered magnetic field 
and contrast it to the previously studied uniform field case. We find that, depending on parameter 
values, a staggered field can provide better conditions for a perfect entanglement transfer, while 
even a modest amount of exchange anisotropy appears to have a strong detrimental effect. We also 
study interactions between different waves of entanglement and assess the possibility of simultaneous 
transmission of multiple bits of quantum information. 



A successful implementation of quantum computing 
protocols hinges on the possibility of achieving a high- 
fidelity transfer of quantum states between the different 
parts, such as 'core processor', 'storage', etc., of a quan- 
tum register. 

To that end, a recently proposed variety of 'all-spin' 
architectures with 'always on' interactions has the po- 
tential of circumventing both the problem of conversion 
between physically different degrees of freedom and that 
of controling inter-qubit couplings individually. 

Specifically, the recent proposals have focused on the 
use of one dimensional chains of permanently coupled 
spin 1/2-like variables ('stationary qubits') whose prac- 
tical implementation might be possible with cold atoms 
in optical lattices, Josephson junction and quantum dot 
arrays, or other two-level systems. 

Being a salient feature of any periodic arrays, the 
translational invariance of spin chains facilitates an emer- 
gence of quantum states with definite momentum ('flying 
qubits'), thereby enabling a high-fidelity spatial transfer 
of entanglement which is viewed as a vital resource for 
quantum computing. 

In most of the proposed designs, the local terms in the 
underlying Hamiltonian ('magnetic field') which can be 
used to perform one-qubit gates are generally expected 
to be much stronger than non-local 'exchange couplings' 
implementing two-qubit ones. Studied thus far are spin 
chains with ferromagnetic 1 , antiferromagnetic 2 , short- 
range and long-range, as well as SU(2)- versus U(l)- 
invariant exchange couplings. On the other hand, much 
less attention has been paid to different patterns of the 
local field, almost all the studies being conducted in the 
(from the technical standpoint, simplest) case of a spa- 
tially uniform field. 

Considering, however, that a local field can be used for 
both, initialization and readout of a quantum register, a 
comprehensive analysis of different local field distribu- 
tions is strongly warranted. To that end, in this paper 
we study the time evolution of entanglement in the pres- 
ence of a staggered field and contrast it to the previously 
explored uniform field case^. 

As an example of the system for which quantitative 
predictions can be made, we use the exactly soluble gen- 
eralized anisotropic XY-model described by the iV-spin 



Hamiltonian 



N 



H = Xj2[a+l)SfSf +1 + a-j)SVSV +1 }-J2biS!. (1) 

i=l i=l 

where 7 is the anisotropy parameter ranging between 
between (isotropic XAT-model) and 1 (Ising model). 
Furthermore, A = J/B is the nearest-neighbor exchange 
coupling and 6; = Bi/ B is the normalized local magnetic 
field in the z-direction, both in units of the average on- 
site field B. In what follows, we focus on the staggered, 
bi = (— l) 1 , field configuration and contrast it to the pre- 
viously studied uniform, 6» = 1, one. 

The overall sign of the exchange term in Eq.(l) is cho- 
sen to be positive, although it can be readily inverted 
under a 7r-rotation Z = Ili ^2i+i( 7r ) a bout the i-axis 
in the spin space on every other site: S^l^L 
ay . qv q x >v o x ,y a x ,y , q x -y cz . oz 

°i ~ * i ' °2i+l ~ * °2i+l> D 2i ~ * °2i > °i ~~ * °i 

Likewise, under a 7r-rotation X — EL -^fi+i M about 
the i-axis the Hamiltonian of the XY-mode\ in a stag- 
gered transverse field transforms into that in a uniform 
one, albeit at the expense of the signs' of the x- and y- 
spin couplings becoming opposite (note that in the Ising 
limit the two models become identical). Moreover, at ar- 
bitrary 7/O the rotated Hamiltonian can be brought 
back into the form of that in a uniform field by a simple 
rescaling: XH (A, 7, (-l) 1 )^ 1 = H(Xj, 1/7, 1). 

Instead of resorting to the above (ill-defined in the 
isotropic limit 7 = 0) identity, we diagonalize the rotated 
Hamiltonian directly by applying the Jordan- Wigner 
fermion representation of the spin operators^ 

at = 2c\ Cl -l,af = Htf)ci, a- = [1(^)4- (2) 



j<i 



followed by the Bogoliubov transformation in the mo- 
mentum space r/k = UfcCfe — ivf~c_ k , with u 2 k + v\ = 1, 
as a result of which the Hamiltonian takes the form 

H = Ek^kivlvk - |). 

In terms of the functions A). = 1 + A7 cos(fc) and B^ = 
Asin(fc), the Bogoliubov coefficients and fermion disper- 



l-A k /u> k 



, and 



sion relation read - 

In the position space, the time evolution of the fermion 
operators is given by the equation of motion dej (t) / dt = 



2 
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FIG. 1: Concurrence and fidelity , i sites re- 

moved from the initially entangled site at (0, 2) for the stag- 
gered field model with 7 = 0, A = 0.05. The dashed lines 
represent the results obtained for a uniform field. 



FIG. 2: Concurrence d t i+2 and fidelity -Fi.i+2 i sites away 
from the initially entangled site at (0, 2) for the staggered 
field model with 7 = 0, A = 0.5. 



i[H,Cj(t)], whose solution is Cj{t) = J2i a i-j(t) c i 
i/3i—j(t)c, with the time-dependent coefficients 

a j( f ) = wJ2 k cos(kj)[cos(u k t) - sw(w k t)], 



(3) 



For 7 = and A< 1 these expressions assume the form 



a 3 {t) 







j even 
j odd 



(4) 



/??'(*) ~ W Efc sin(fcj) sin(fc) sm(ui k t) 



The above formulas allow one to determine the time 
evolution of an arbitrary initial state, although if the ini- 
tial state is highly mixed this process becomes computa- 
tionally rather cumbersome. Therefore, in line with the 
previous research^ we focus upon simple initial spin con- 
figurations where all of the spin states are aligned with 
the field, with the exception that some pairs of spins be- 
gin in an entangled state. 

In the staggered field model, we first look into the be- 
havior of a state with only one singlet 



hi 



(4 + c,)]luu 



(5) 



formed between the spins i and j residing on the same 
sublattice. 

The entanglement properties of this state are described 

(2) 

by a two-spin reduced density matrix p m \n which can be 
obtained from the density matrix of the entire system 
with all the sites but i and j traced out. In the basis 
{|TT> , 144) , lit) , 144)} it can be cast in the form: 



.,(2) 



(*) = 



/ (p+p+) (p+ a ~ 



(P^PZ) 

(p^i) 



\(7 ra P n ( 

ip-p- 

\ m n 



(6) 



where P±(t) = ±(l± <(*)). _ 

Since tpij and H are both invariant under parity, the 
expectation values (<r^), (of„), (cr^O, and 
(«) vanish. 

The remaining, non-vanishing, equal-time correlation 
functions appearing in Eq.(6) are most easily evalu- 
ated by rotating with the operator X, so that the field 
becomes uniform and the initial state takes the form 
^,i = ^(cl(0)-c](0))|^). 

Thus it is necessary to evaluate matrix elements of the 
form (^iXFiP^P+^^a^X- 1 ]^) which can be 
computed directly via Wick's theorem, or by a method 
involving sums of Pfaffians^. 

As a convenient measure of pair-wise entanglement we 
will use concurrence defined according to the formula 5 

dj = Max {0, A 4 - A 3 - A 2 - Ai} , (7) 
(2) (2) 

where Ai are the eigenvalues p\ ® o- y {p\ j)*<J y ® a y in 
ascending order. For a comparison, we will also compute 
fidelity defined as follows: Ff - = Tr{p^J \tp) (tp\), where 
\ijj) is the singlet Bell state on the initially entangled pair 
of spins. 

Starting out in the initial state (5), in the isotropic 
(7 = 0) strong field (A <C 1) limit the system under- 
goes a time evolution described by Eq.(4) which can be 
viewed as the propagation of a singlet (EPR) pair which, 
after having been initially created at the sites i and j, 
starts moving in both directions away from these sites at 
a speed A 2 /4. 

Thus, consistent with the general analysis of Refs^, 
the entanglement propagation can be characterized in 
terms of the dispersion and damping (apparent decoher- 
ence rate) of underlying elementary excitations, as long 
as only a limited number of which is generated in the 
process of a quantum state transfer. 

In the position space, the operator of a wave packet 
corresponding to each of the two spreading entanglement 



3 



waves is a superposition 

djK^iV-Mj^wiXH/W- (8) 
I 

which receives its largest contribution from the sites on 
the same sublattice as i, thereby manifesting the fact 
that the entanglement wave remains largely confined to 
one sublattice. 

Figs. 1 and 2 show the results of the calculation of 
and F^i + 2 for the isotropic (XX) staggered field 
model with A = 0.05 and A = 0.5, respectively. In the 
plots, the initially entangled pair is on the sites num- 
bered and 2, labeled as (0,2), while the fidelity and 
concurrence are computed on the sites (2,4) and (6,8). 

In the strong field limit, the quality of entanglement 
transport in the staggered field model approaches that 
of the uniform field one (the dashed lines). The appar- 
ent 'ripples' in the plots of Fig.l should be attributed to 
vacuum fluctuations arising due to the fact that for any 
finite A the fully polarized (classical Neel) initial state 
differs from the true ground one. 

With increasing A the vacuum fluctuations become 
more prominent and, concomitantly, the amount of trans- 
mitted entanglement decreases, as seen in Fig. 2. While 
the singlet fidelity remains approximately constant, this 
reduction of the entanglement transfer is due to the emer- 
gence of additional states of the form a \H)+b |j.T), whose 
contribution is governed by the coefficients /3;_iA from 
Eq.(4). 

Next, we study the time evolution of an initial state 
with two individual entangled EPR pairs 

= ^(cl+c i -(c]+ Cj -))(4+c fe -(cj+cO) inn ...> (9) 

Specifically, in Fig. 3 we contrast the two cases when both 
pairs reside on the different (top picture, (i,j),(k,l) = 
(24, 26), (47, 49) and same (bottom picture, (i,j),(k,l) = 
(24, 26), (46, 48) sublattices, the brighter areas indicating 
greater values of the nearest sublattice-neighbor concur- 
rence. 

The outward-moving entanglement waves propagate 
away from the initially entangled sites with the same 
speed A 2 /4 as that in the one-singlet case. The lack of any 
detectable concurrence among nearest lattice-neighbor 
sites in the sublattice not containing the initially entan- 
gled sites is a direct consequence of Eq.(4), in the strong 
field limit. 

The other two inward-moving wavefronts collide and 
then either pass through or annihilate each other, de- 
pending on which sublattice they reside. If the initially 
entangled sites are on different sublattices, then accord- 
ing to Eq.(4) all of the fermion operators Cj(t) that con- 
tribute to the time evolution of the system will have 
odd indices for one wave and even indices for the other. 
Therefore, in this case the time evolution of the two col- 
liding waves proceeds independently. 

If on the other hand, the two waves are on the same 
sublattice, Eqs.(4) and (8) may be used to show that, 
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FIG. 3: Time evolution of C M+2 for A = 0.05 and 7 = 0. 
The 2 initially entangled pairs of sites in the top plot are on 
different sublattices. In the bottom plot the two pairs are on 
the same sublattice. 



when the waves intersect, all entanglement they were 
transmitting is exactly negated by the other wave. This 
negative interference is a generic property of the uniform 
field models where it can present a serious problem for a 
simultaneous unimpeded propagation of multiple flying 
qubits. 

As A is increased the term in Eq.(4) proportional to the 
function (3j (i) not only serves to dampen the transmission 
as already mentioned, but also to diminish the degree to 
which the waves confined to different sublattices maintain 
their integrity. Formally, this occurs due to the function 
Pj(t) taking a non-zero value for all values of j, which 
allows for the two waves to mix, resulting in their loss of 
coherence. 

The observed nearly exact cancellation of the trans- 
mitted concurrence is not a general result, since it only 
occurs when the two colliding waves are created at pre- 
cisely same time on the same sublattice. However, the 
ability of entanglement waves to pass through each other 
if they propagate on different sublattices is independent 
of the strength of the entanglement carried in each wave 
and also of the actual state that is transported. 

Finally, we investigate the effect of a small anisotropy 
(7 <C 1) in the strong (staggered) field limit (A < 1) 
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FIG. 4: Time evolution of the concurrence Ci,i+2 for the 
anisotropic XY system in a staggered field, 7 = 0.0125 and 
A = 0.05. 



where the fermion dispersion becomes uik ~ 1 + + 
7Acos(fc) — ^-cos(fc). The mixed cosine terms do not 
allow for a closed form for the function aj (t) unless A 3> 
7, in which case Eq.(4) still holds. 

In the opposite limit (A <C 7), one finds that (3j(t) is 
again negligible, while the result 



i(t) w i~ j Jj(j\t). 



(10) 



suggests that the initial excitation is no longer confined 
to its sublattice but propagates over the entire lattice, the 
corresponding wavefront traveling with the speed 7A. 

The intermediate regime A ~ 7 <C 1 is presented in 
Fig. 4 where the nearest sublattice-neighbor concurrence 
is plotted for only one side of the propagating wave in 
the case of 7 — 0.0125 and A = 0.05. The plot shows 
a presence of two types of elementary excitations, one 
propagating on all sites at a speed 7 A and the other, 
with a speed A 2 /4, which is confined to one sublattice, 
as in the isotropic limit. 

As the degree of anisotropy increases, there is a trans- 
fer of spectral weight from the latter to the former mode, 
which starts to dominate the transport of entanglement 
already at 7 ~ A and continues all the way to the Ising 
limit. 



In conclusion, we studied the dynamics of the spin 1/2 
XY model in a staggered transverse field for a few simple 
out of equilibrium states. This previously unexplored 
problem is of interest in those prospective designs of a 
quantum register where the local field can be created by 
a ID antiferromagnetic crystal (which systems happen to 
be more abundant than their ferromagnetic counterparts) 
placed in the vicinity of the II? qubit array. 

The use of non-uniform (especially, alternating) fields 
in qubit chains might also be advantageous in light of the 
their ability to create additional entanglement^. More- 
over, the study of various field configurations is impor- 
tant for optimizing local field profiles^ as a part of the 
general problem of finding the best parameter values for 
artificially engineered qubit chains (most notably, cold 
atoms in optical lattices or Josephson junctions arrays). 

To that end, we found that in the strong field limit the 
system approaches the uniform field model in its abil- 
ity to transport entanglement, however, the speed of the 
state transfer appears to be lower because of an extra 
factor of the inverse field strength, which factor needs to 
be small in order to transport states with an acceptable 
degree of fidelity. 

In this regime, the two sublattices of the chain act inde- 
pendently, so that waves traveling on one sublattice will 
not interfere with those on the other one. After their 
collisons, the integrity of such entanglement waves re- 
mains largely intact, thereby facilitating the possibility 
of transporting multiple quantum bits at the same time. 

On the other hand, a generally slow speed of transmis- 
sion in the staggered field case increases the propagating 
state's exposure to any environmental sources of deco- 
herence, thus manifesting yet another potential tradeoff 
which needs to be accounted for in the prospective de- 
signs of a practical ID quantum register (of course, the 
emergence of mutually contradicting criteria is a rather 
common place in the field quantum computing in general 
and in the qubit-chain designs in particular 9 ). 

Lastly, by examining the effect of anisotropy in the 
XY staggered field model, we find that the transmission 
starts to lose the characteristics of staggered field prop- 
agation already for moderate values of the anisotropy 
parameter 7, as compared to the ratio of the coupling 
strength to the field strength A. 

Our results provide a further insight into the non- 
equilibrium dynamics of qubit chains and facilitate an 
ongoing quest for the optimal conditions of a reliable 
quantum information transfer. 
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